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I . OBJECTIVE 



In this thesis the Sequential Rayleigh Test and three 
sequential binomial tests, which are applicable to testing 
bombing system accuracy, are compared by computer simulation. 

The objective of this thesis is to investigate expected 
sample size, variance of sample size, error rates of these 
tests, and to investigate their robustness. 
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II. DESCRIPTION OF THE TESTS 



A. INTRODUCTION 

There are two types of tests of a system (which is 
taken to be a bombing system in what follows) . The sample 
size is determined as a direct result of the experiment 
in one case (and is therefore random) , and the sample size 
is selected prior to commencing the test in the other case. 

The former procedure is called a sequential test [1] , 
The test procedures being considered here are sequential, 
and are based on assumed circular normal distribution (that 
'is, bivariate normal) with the same variance in each co- 
ordinate. In the coordinate system of the target plane it is 
assumed that: 

X-^ Normal ( 0, a^) 

— ■ Normal ( 0, a^) . 

The origin of coordinate system is the target, and x arid 
Y are distances from the weapon impact to the target along 
the X and Y axes . 

If it is further assumed that X and Y are independent, 
then — ^ is normal with mean 0, variance 1. -~r- is normal 

with mean 0, variance 1. Thus ( ^ has Chi- 

square distribution with two degrees of freedom. 
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The density function of Chi-square with two degrees of 
freedom is [2] . 



(t) 



1 

2-r(i) 




t>o 




-t/2 



This is the exponential density function with parameter 
X— ^ . Here/ 

X^+Y^ , where X^ + Y^ is squared miss 

2 2 ~ 2 
o o 

distance of the impact from the target. Let this be another 
random variable Z. Then, 

Z ^ Exp (i) 
a2 

The density function of Z is derived as follows: 

F^(z) = P 

= P 



[ Z< z ] 




a o 



8 



Let 



be Z ; 



then FgCz) = { -^) 



Where Z ^ is exponentially distributed with parameter \= ^ . 



Then 



F = 1-exp [-Z/2. cr^] =F (z) . 



From the relation that — =7,' the density function 

cr"^ 



of Z is [5] 



2a 



exp[-z/2.a^] = 



which is the exponential density function with parameter X= 1 



2a‘ 



Assume C is the median of Z (C represents CEP which is de- 
fined to be the median of the squared Circular Error Probable 
[ 6 ] ) . 

Then by definition: 

F„(C) = i 



1 -exp[-C /2.a^]= ^ 



C 

2a' 






a 



C 



Thus Z ^ Exp ) 
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If a bombing system has been specified to have median radial 

miss distance Yo and if a system with radial miss distance 

2 

Yi is unacceptable, this can be tested with Co= Yo as the 

2 

null hypothesized median and Ci =Yi ss the alternative 
hypothesized median. Here, in2 ln2 are the parameters 

Co ’ 

under null and alternative exponential distributions, 
respectively . 

In what follows Co will be treated as median under null 
hypothesis and as median under alternative, where these 
quantities relate to the squared radial miss distribution. 



B. TEST CASE A; SEQUENTIAL RAYLEIGH TEST 

If a system is to be tested with null hypothesis 
H^; CEP ^ ~ alternative hypothesis H^^; CEP ^ ~ 

with type I and type II error rate a and g, then a sequential 
test for H^ against H^^ can be defined as follows. We shall 
call this the "Sequential Rayleigh Test" in what follows; 
it is an application of Wald's Sequential probability ratio 
test to the exponential situation described above. 

Two positive constants are chosen, B and A, where B < A. 

At each stage (n-th experiment or observation) , the probability 
ratio n f(Zj, Cj^) ^ called the likelihood 



n 

i=l 



f(Zi,* Cq) 



ratio [3], is computed, where f(z^^; Cj) is the exponential 

density function with parameter £ n2 , j = 1.2 

C. 

D 
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<A, then another observation is made (this 



n f(z., Ct) 

ill £(zi, c"T 



St 



means the test enters the (n+1) stage) . 

If the likelihood ratio does not fall in the interval (B,A) , 
called the continuation region, the test terminates. 

In termination, the conclusion is to; 



n 



f(z^, 



Accept if C T ) ~' i 



n 

Reject if II 



f(Zi/ 
f(z., C ) 



> A. 



i=l ’ i^ o’ 

For a test with approximate level of significance a and power 
1-3, one may define [1] 



A = ' and B = 

a 



3 

.-a 



Thus, an approximate bound for each stage can be obtained as 
follows : 

The explicit form of the likelihood ratio is 



n f(z.,C,) 
i"i 

_ {5-n2/C-|^) .exp [-z^.Hn 2 /C 2 ^] . - - - . (I!,n 2 /C 2 ^) . exp [-z^. J2,n2/Cj^] 
{in2/c^) .exp [-Zj^ . Hn2/c^] . - - - . (5,n2/C^) .exp[-z^.£n2/C^] 

(IIn2/C,)^. { p z.) .exp[-£n2/C,] 

= ^ i=l ^ ^ 

(£n2/C)^.( Z z.) .exp[-£n2/C ] 

o i^l 1 o 
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This is equated to A and the logarithm is taken to obtain 



rejection bound for the n-th stage; 



i^n(-i^)=n.tn/^ 

1^1 



n 

I Z..1 
i=l 



“1 



ln2 - Jin 2 



Cl 

1 o j 



The test rejects if 



E z,> -dn/l-eN 'l = R„ 

i=i 1-C77 " 



n 



where E z- is sum of squared radial miss distances, 
i=l 



The same procedure can be applied to B; 




The test accepts 
n 

I 2 
i=l 

(Jln2) . 





= \ 



Here, it is assumed > C^, which is, we envision, 
in the bombing system test. 



true 



C. CASE B; SEQUENTIAL BINOMIAL TEST WITH Pq = 7 

In the study of bombing system the target is defined to 
be a point on the impact plane and impact of a bomb within 
(over) some distance /r from the target is defined to be a 
hit (miss) . 

2 

A null hypothesis that CEP = is to be tested against 

2 

an alternative hypothesis that CEP = with Type I and 

Type II error rates a and B, where 

2 

Also this system can be tested with CEP < as null 
2 

hypothesis and CEP ^ as alternative hypothesis without al- 
tering the test procedure. 

Let Pq be defined as the probability of hit under the null 

hypothesis, and to be the hit probability under the alterna- 

2 

tive hypothesis. Then under CEP < C^, which says the 

true median of squared radial miss dis tence is less than or 
equal to C^, implies the probability of hit is greater than 
or equal to 0.5. 

Similarly, CEP > implies the probability of hit is 
less than or equal to 0.5 under the alternative hypothesis. So 
the new hypotheses are defined as P^> 0.5 and P^^ < 0.5. 

From this hit or miss criterion r, a value of squared radial 
miss distance may be obtained which gives P^ value of 0.5; 
i.e. , if 
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-Zn2 . r 



F^(r) 



1-e 




0.5 



then S,n(0.5) = (Jln2) . (- — ) t so r = C^. 

The alternate hit probability is found as follows: 

?! = [hit I CEP^ = C^] 

= P^ [observed squared radial miss distance < r | CEP -C^^] 




= 1 - e 



ln2 

^1 



= 1-2 



Let the random variable be defined as: 



z . = 
1 



0 if miss (squared miss distance > 

1 if hit (squared miss distance < 
Then the likelihood ratio becomes 



n 

n 

i=l 



TTITTp^ 






( 1 / 2 ) ‘ ( 1 / 2 ) 



n-Z 



^1 

.1^ 



(1-Pl) 



n 



,n 



r) 

r) 



n 

Where Z = z z- 
i=l ^ 
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For a sequential probability ratio test for the binomial 
situation, the experiment is continued as long as this value 
remains between B and A. 

An approximate acceptance boundary is found by substituting 



g for B [1] . 
1-a 



Then Jin 



■g 



n 






i=l 



Substituting (1-2 *“ 0 ^*^!) for gives 



Jin 






• In 






n 



Solving for z^, which is a convenient test statistic; 



n 

Z 

i=l 




Jin /-g 






-n . ( 1 - Cq/C^ ) • Jln2 



in [2 ° ^ - 1) 



A 



n 



Cq/Ci 

(The inequality changes because implies 2 - 1 < 0.) 

Where ^ 2 ^ is the number of bombs which hit the target, out 

i=l i_g 

of the total fire, n. And substituting — ^ for A, an 

approximate rejection boundary is found to be: 



15 




The test now operates as follows: 



in stage n. 

Accept Hq if j ^i - ^n' 

Reject Hq if j R^; 

i=l 

Continue to stage n+1 otherwise. 



D. CASE C: A SEQUENTIAL BINOMIAL TEST WITH NULL PARAMETER 

WHICH MINIMIZES E [N] UNDER Hq. 

2 

This case also has null hypothesis CEP <Cq and alter- 
2 

native CEP >Cj^. Let PQ(r) be the probability of hitting a 
target of radius \Tr^ and let P^Cr) denote that probability 
under H^^. In mathematical form: 

j!,n2 . r 

2 ^ n 

PQ(r) = PQ(Z<r|CEP = Cq.) = 1-e 



Pl(r) = Pj.(Z<rlCEP^ = C^) 



1-e 



Jln2 




r 
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c /c 

Thus: Pj^(r) = 1 - [1 - PQ(r)] ® ^ 

It is desired to determine r so as to minimize expected sample 
size n, required to test Hq vs with Type I, II error rates 
a, 3 respectively, using the Binomial Sequential Probability 
Ratio test. 

The average sample size function is 



^ ... (l-L(P) ) . JlnA + L (P).JlnB 

^ 



ti] , 



P.£n 






+ (1-PUn 




[4] 



where P is the true probability of hit and L (P) =P^ [accept Hq|P]. 
If PqCt) is chosen to be the underlying hit probability, then 
LCPpCr) )=l-a. 

By substituting L (Pq ( r) ) =l-a , 

A = 1-3 , P3^(r)=l-(l-PQ(r) )-Cq/Cj^,B = 3 

a 1-a 

in E[N] function, the expression above becomes 



E(N) 



aJln 



(k±] IJA 

\ a. / + (1-g) £n(,l-g/ 



l-d-P.(r)) ^ -l+C./C, 

' ] + (l-PQ(r))Jln[a-Po{r)) 



Pq (r)j[Jln 



Po(r) 



If a value of PQ(r) which minimizes E(N) is obtained, r is also 
obtained from PQ(r). The numerator is a negative constant as 
long as 3<0.5 and a<0.5 and the denominator is a negative variable 
depending on r. 
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Therefore, in order to minimize E(N) it is necessary only to 
maximize the absolute value of the denominator. 



Let 




_1 

k 



Then for various values of K, ^^ir) may be found. 

Figure II- Is hows a plot of PQ(r) vs k resulting from this 
minimization. As k approaches infinity PQ(r) approaches 1.0; as 
k decreases to zero PQ(r) decreases to around 0.63. Specifically, 



for k = 2 

PQ(r) = 0.8416, P^(r) = 0.6020. 
From this r is found to be 

- = 2.65836.C^ 

Let Z be a random variable such that 
2 _ if squared miss distance 

(i if squared miss distance 



> r 
< r 



where r = 2.65836. 

Then the hypotheses in this Binomial test becomes: 



Pq = 0.8416, : P^^ = 0.602. 

If Pq = 0.8416, Pj^ = 0.602 are substituted in the likelihood 
ratio (similar to Case B above) , acceptance and rejection bounds 
at stage n are; 

Accept if 



ill 2^ > -1.2564 



• rr77) + n. 0.7333 = A^ 

-L vX 11 
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FIGURE II-l 
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I 




Reject if 

? Z. < - 1.2564.J?n ( + n. 0.7333 = R„ 

. , 1 a n 

1=1 

Continue otherwise. 



E. CASE D; A Sequential Binomial test with null parameter 

which minimizes maximum value of E(N) 

In the binomial situation described above, the hypotheses 

about the parameter are 

H_; P = P 
o o 



P = 



1 * 



The maximum average value of E(N) occurs very nearly at 



I n 



P- = 



I n 



l-Po(r) 

l-Pl(r) 



Pl(r) 



in 



l-Pl(r) 



[4) 



from which 
E(N) = 



~ 5AA. £ n B 



SA1 



'Pl(r) 



— £4 



'l-P^(r) 

T^P^TFT 



ihA. inB 





l-d-P^(r) fo^^l 


I 

. iln 


a-P^(r) f o/ci ' 


' 




Po(r) 

1 ^ 




l-Po(r) J 





(P^(r) = 1 - [1-P^(r)f q/^ 1, See Case B) 

The numerator is a negative constant for given a and $. 
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Similarly as in Case C by maximizing the absolute value 
of the denominator/ the maximum value of E[N] is minimized (or 



at least nearly so) . 

Figure II "3s hows the relation between and k = 




Specifically for k = 2 

PQ(r) = 0.89867, P^(r) = 0.68167 with r = - Z n(i-Q. 89867) = 



3.30286 is obtained. 

Hence the hypotheses for this case are P=0. 89867, 

P=0. 68167, and "hit" is defined by squared radial miss 
distance less than or equal to 3.30286. 

Let Z be a random variable such that 
('O if squared miss distance > 3.3 0286 

Zi = 

II if squared miss distance < 3.30286. 

Proceeding similarly as in Case B, the decision at stage n is 

Reject H if 
o 



n 




+ n. 0.80552 = R. 



n 



Accept H if 



o 



n 



E Z. > -1.42167 . An 
i=l ^ 



6 ) + n. 0.80552 = A. 



a 



n 



Continue to state n+1 otherwise 
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I 




III. DETERMINATION OF SIMULATION FACTORS 



For convenience the hypotheses 





CEP^ = 1 




Hi: 


CEP^ = 2 


were selected. 


This means k = 2 


in Case C and D. 



a 



Type I and Type II error rates were selected to be 0.05. 

But actual simulation gives Type I error rates .0258, .0468, 

.0396, .0292 and Type II error rates .0396, .0402, .0420, 

•0400 for Case A, B, C, and D, respectively. This point is ex- 
plained in Chapter 3.3 of [1]. Approximate error rates of 

0.05, 0.05 (Type I, II) were obtained by adjusting the bounds 

1 — 6 

A and B, which is possible by changing a, 6 in A = — ^ , B = 

The adjustment factors used here are: 

CASE A a= 2. Ox 0.05 $ =1.2x0.05 

CASE Ba=1.2x0.05 6 =1.2x0.05 

CASE C a= 1.5 X 0.05 B = 1.15 x 0.05 

CASE Da=1.5x0.05 3 =1.1 x0.05 

How many replications are enough? Assuming N^, N^ have Bi- 
nomial distribution with probability of success 0.05, a sample 

A. /\ 

size n is found such that = 1- a , where a is the 

significance level, P^, i=l, 2 is the estimation of P^^ . The above 
equation implies 

/V. /N 

Pj, (-C < “ ^2 < c) = 1-a 

Pr(-c < 5^1 - ^^2 < c) = 1-a 
n n 

P (-nc< N - N <nc) = 1- a 
r 12 
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By the Normal approximation , - N 2 )-~ Normal with mean o and 

Variance 2np(l-p). 



Thus 



[ - 



nc 



\^/2np(l-p) 



% - ^2 
V'2np (1-p) 



nc 



v' 2np (1-p) 



] = 1-a 



gives 



nc 

2np(l-p) 




a 

7 ~ 



The following table shows n for various values of a and c. 





0.2 


0.1 


0.05 


0.025 


0.02 


c Y 


1.285 


1.645 


1.96 


2.24 


2.33 


0.2 


3.9 


6.4 


9.1 


11.9 


12.8 


0.1 


15.7 


25.7 


36.1 


47.6 


51.5 


in 

o 

• 

o 


62.7 


102.8 


14 5.9 


190.7 


206.3 


0.025 


250.9 


411.3 


583.9 


762.7 


825.2 


0.01 


1568.7 


2570.7 


3649.5 


4766.7 


5157.5 


0.001 


156866.4 


257072.3 


361237.5 


476672.0 


515745 .5 
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Assuming c = 0.01 and a= 0.02, an approximate value o£ 
5000 is obtained through the table. By this number of re- 
plications, obtaining a difference in estimated error rates 
greater than 0.01 is significant at level 0.02. Exponential 
random samples were generated by the Monte Carlo method. For 
an exponential variate T to have median m, it is necessary 
to use the scale parameter \ = I n 2 

n O ^ 

- t 

SO F.p(t) = 1-e 

But U = F-p (T) is uniformly distributed on (0,1) [8]. 

-T 

Thus 2 = l - f(T)= 1-U is also uniformly distributed. 

Finally t ,. Zn2 = -£n U, or 
m 

T = - m. 2, n U 
in2 

By changing the median of the population sampled, the 
operating characteristic function, expected sample size at 
termination, and its variance can be estimated for each test, 
and those can be compared. 

The medians to be generated are those which result in a 
0.1 difference of operating characteristic function values 
in Case C, which is based on minimizing the maximum E(N), 
those which yield Type I, II error rates, three points which 
yield approximate maximum E(N) value (for Cases B, C, D) , and 
five more points in both tails. 
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By changing the skewness of the sample distribution, 
comparison of the robustness of the four tests of hypotheses 
about the medians under null and alternative hypothesis was 
performed. 

This is based on the assumption that the underlying 
distribution is WEIBULL with shape parameter a and scale 
parameter X . 

If T ^ WEIBULL{a, X }, 



the distribution function of T 



F^(T) = 1 - e 



-(XT)“ 



is 



Let m be the median. Then 



a 



Pj. [T > m] = [T < m] ; 

1 - 0.5, 

Jin(o.5) = - (X.m) 

Solving for X , 

X = 

Now 1-e . — - u(o,l), so 

£n(l-u) = - (XT)“ , 

T = y [Jln(U) ]'^ , where U 

Substituting X = i 

m 

1 , , . 
T = m . .jm,u.''^-WEIBULL{a (m2 



U (0,1) . 
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Following histograms in Figure III-l show the effect of 
changing a in WEIBULL. As a decreases the distribution is 
widely spread and it is said that the distribution has heavy 
tail (Figure 111-1/ ®). In the opposite case; i.e. a increases, 
it has light tail (Figure III-l, A) . 

Figure III-2 is observed keeping a fixed at 1 and median is 
0.5, 1.0, 2.0 in WEIBULL. 

This is exponential distribution with median 0.5, 1.0, 2.0. 
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FIGURE III-2C, WEIBULLU, 0.3466} 



IV. RESULTS AND COMMENTS 



Simulation of the four tests were performed 5,000 times 

2 

(5,000 replications). In each replication CEP =1; 

2 

CEP =2; and Type I, II error rates are 0.05. 

The following tables IV-1,2,3 show the mean number of 
terminations, their variances, and number of acceptances of 
Hq in 5,000 replications. 

Table IV-lshows the effects of changing median values, where 
the shape parameter of the WEIBULL variate was fixed at 1.0 
(the exponential variate with changing median) . Table' IV-2 and 3 
shows the effects of changing shape parameter, the median being 
fixed at 1.0 in table IV- 2and at 2.0 in table IV-3. 

For ease of comparison these were graphed. Each figure 
contains four curves and each curve represents the case A, B, 

C or D. 

In testing the difference of E [N] , null hypothesis would be 

E[Nj^] = E[N 2 l/ where Nj^, N 2 are random variables representing 

the number of termination by test cases which one wishes to 

test the difference of the E[N]. 

By the central limit theorem [2] the difference in E[N] 

greater than 1.96 x d is significant at level 0.05. d is 

n 

estimated standard deviation of N. But it is different 
between the test cases being considered. 

The larger one might be selected so that we may not 
reject erroneously. 
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THIS CASE IS C6SERVED CHANGING MEDIAN ANC SHAPE PARAMETER IS FIXED AT 
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TABLE IV- 



lb UtibtHVtU LHANOlNCi bHAPfc HAKAMETER ANC f^tOlAN IS FIXED AT 
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REFERENCE FOR THE FIGURES IV-1 - IV- 9 



CASE A 



+ + + CASE B 

* * * CASE C 

X X X CASE D 



FIGURE 
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Expected 
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FIGURE 


CN 
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Variances 
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FIGURE 


IV-3 , 
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Operating 


Characteristic Curves 
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FIGURE IV-1 
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The interpretation of the figures is as follows; 

a. Variances or standard deviations are approximately pro- 
portional to expected sample sizes. 

b. Test Case A contains the minimum expected sample size and 
Test Case B includes the maxim\im expected sample size in 
almost all instances. 

c. Case C and Case D yield similar results when median is 
changing and shape parameter is fixed at 1.0. 

d. Test Case B is invariant about the change of skewness 
when the underlying median is 1.0 (See Figures IV-4,5,6). 

e. Case B also gives the lowest error rates if the bomb im- 
pact has a heavy tail distribution under the alternative 
hypothesized median, and has light tail distribution under 
null hypothesized median (Figure IV-6.9). 

f. Case C gives the lowest error rate when the true median is 
2.0 and bomb impacts are clustered around the median. 

NOTES ; 

a. The adjusting factors of the error rates table were devel- 
oped iteratively and there is no guarantee that these ad- 
justing factors are the best ones. 

b. There could be further investigation of effects under er- 
ror rates other than 0.05, or of changing null and alterna- 
tive hypothesized median values. 

c. The Fortran program for this simulation is attached as an 
Append ix . 



47 



r 

i 




V. CONCLUSION 



If it is certain that the bomb impact is Rayleigh distributed, 
then Sequential Rayleigh Test is appropriate to test the system. 

In case one has doubt about the bomb impact distribution, 
the Sequential Binomial Test gives better tests than the Se- 
quential Rayleigh Test. 

Further, if there is some reason to believe that the 
median is likely to be 1 (so the test is likely to accept H^) , 

Case B (Null parameter P^=0.5) is better than any of the other 
three tests. 

If the median is not likely to be 1, test case A (Sequential 
Rayleigh Test) or test case C (which has null parameter that 
minimized E[N] under H^) are the best tests. 
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APPENDIX 



DIMENSION BCUN0{300t2»4),0ATA(4»4) ,T EST { 5000 1 2 , 4 ) 

BOUND CONTAINS DECISION BOUNDARIES FOR EACH STAGES 

TEST FIRST ELEMENT REPRESENTS NUMBER OF REPLICATION 

SECOND ELEMENT REPRESENTS DECISION 

THIRD ELEMENT REPRESENTS TEST CASE(A,B,C OR D) 

DATA TEST RESULT, FIRST ROW MEAN NUMBER OF TERMINA 

TION SECOND ROW VARI ANCE, THIRD ROW NUMBER ACCEP 

TANCE FOURTH ROW NUMBER REJECTION 

ISEED= 688777 
NC0UNT=50C0 

READ (5, 2) ALPHA, BETA, CO, Cl 

CO MEDIAN UNDER THE NULL HYPOTHESIS 

Cl MEDIAN UNDER THE ALTERNATIVE HYPOTHESIS 

WRITE! 6, 140) ALPHA, BETA, CO, Cl 

140 FORMAT(5X, *ALPHA=' ,F6.3, • BE TA= • , F6. 3 , • C0=',F6.3 

’i',* C1=',F6.3) 

R=-C0=*AL0G( .1584)/( AL0G(2. ) ) 

R HIT MISS CRITERION IN CASE C 

R2 HIT MISS CRITERION IN CASE C 

P40=.8S8657 

R2=( ALOG( 1-P40) )/AL0G(.5) 

P41=l.-l./2*=<‘(R2/2. ) 

2 FORMAT! 4F6.3) 

COMPUTE BOUND OF REJECTION AND ACCEPTANCE 
B1=AL0G! 1.2^BETA/!1-2.*ALPHA) ) 

A1=AL0G! ! 1-1.2*BETA) /!2.*ALPHA) ) 
B2=AL0G!1.2-BETA/!l-1.2=i= ALPHA) ) 

A2=AL0G! !l.-1.2=i‘BFTA )/!l .2’i‘ALPHA) ) 

63=ALOG!l. 15=»8ETA/! 1-1.5«ALPHA) ) 

A3=AL0G! !1-1.15^BETA)/!1.5^ALPHA) ) 

84=AL0G! 1. 1*BETA/!1-1.5=^ALPHA) ) 

A4=AL06! ! 1-1 .1«BETA)/! 1.5«ALPHA) ) 

AL=!AL0G!2.) )=M 1 ./CO-i./Cl ) 

AL1=AL0G!2-*!C0/C1)-1. ) 

F1=1-SQRT! .1584) 

AL2=ALCG! .!584>>'P1/!.8416«! 1-Pl ) ) ) 

AL3=!AL0G! 2.) )’*'! 1-!C0/CD) 

AL4=AL0G! ! 1-Pl )/ .1584) 

AL5=AL0G! ! 1-P41 ) /I1-P40) ) 

AL6=AL0G!P41^! 1-P40) /! P4 ! 1-P41 ) ) ) 

DO 5 1=1,300 

BOUND !I , 1,1) = ! I=f'ALOG!Cl/CO) + Bl ) /AL 
BOUND! 1,2,1)=! Al-I A ALOG!CO/C1) )/AL 
BOUND!I ,1 ,2)=!A2 -IaAL 3)/AL1 
BOUND! I,2,2)=!B2-I*AL3) /ALl 
BOUND! I ,1 ,3)=! A3 -IaaL4)/AL2 
BOUND! 1,2 ,3)=! B3-I*AL4J /AL2 
BOUND! 1,1 ,4) = !A4-I *AL5) /AL6 
BOUND! I, 2,4)=! B^ -I «AL5) /AL6 

5 CONTINUE 

300 RSAD!5 ,6 ,ENO=1000 ) FMED, SHAPE 

FMED MEDIAN VALUE TO BE GENERATED 

SHAPE SHAPE PARAMETER IN WEIBULL DISTRIBUTION 

6 F0RMATI2F6 .3) 

WRITE! 6,102) FMED, SHAPE 
KCUNT=1. 

CO 7 1=1,4 
00 7 J=l,4 

7 0ATA!I,J)=0. 

Tl(2,3,4) TEST STATISTIC FOR TEST CASE A!B,C,D) 

DCZN1!2,3,4) — TAKES VALUE 1. IF TEST TERM INATED , OTHE 
RWISE 0. 

10 T1=0. 

T2=0. 

T3=0. 

T4=0. 

N=0. 

DCZN1=0 
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DCZN2=0 

DCZN3=0. 

DCZN4=0 

20 CALL GGUB( I SEED» ItU) 

X = FMED*(ALOG(U) /AL0G( .Sn=^‘*( 1 ./SHAPE) 

N=N+1 

IF(DCZNl.£Q.l.) GO TO 40 
T1=T1+X 

IF(T1.GT.B0LND(N, 1,1) .AN0.T1.lt. BCUND(N, 2,1)) GO TO 40 

TEST(KOUNT, 1,1 )=N 

DATA(1,1)=DATA(1,1)+N 

IF(T1.LE.B0UND(N, 1,1) >G0 TO 30 

TEST(K0UNT,2,1)=-1. 

GO TO 31 

30 TEST(KOUNT, 2,1)=1. 

31 DCZN1=1 . 

40 IF(0CZN2.EQ.l. )G0 TO 60 
IF(X.LE.CO) T2=T2+1. 

IF(T2.GT.B0UND(N, 1,2) .AND.T2.lt .BCUND(N, 2, 2) )G0 TO 60 

TESTCKOUNT, 1,2)=N 

C4TA(2,1 ) = DATA( 2, 1 ) + N 

IF(T2.GE.B0UND(N,2,2) )G0 TO 50 

TEST(KOUNT, 2,2)=-l. 

GO TO 51 

50 TEST(K0UNT,2,2)=1. 

51 CCZN2=1. 

60 IF(DCZN3.EQ.l .)G0 TO 80 
IF(X.LE.R)T3=T3+1. 

IF(T3.GT.B0UND(N, 1,3) . AND.T3.lt. BOUNDIN, 2, 3))G0 TO 80 

T£ST( KOUNT, 1,3)=N 

DATA(3, 1) = DATA( 3, 1) + N 

IF(T3.GE.B0UN0(N,2,3) )G0 TO 70 

TEST(K0UNT,2,3)=-1. 

GC TO 71 

70 TEST(K0UNT,2,3)=1 . 

71 DCZN3=1. 

80 IF(DCZN4.EQ .1) GO TO 90 
IF(X.LE.R2) T4=T4+1. 

IF( T4.GT.B0UND( N, 1,4) .AND.T4.lt. BCUNO(N, 2,4 ) ) GO TO 20 

TEST(KOUNT, 1,4)=N 

0ATA(4, 1I=DATA(4, 1 )+N 

IF(T4.GE.B0UND(N, 2,4) ) GO TO 78 

TEST( KQUNT,2,4)=-1. 

GO TO 79 

78 TEST(KOUNT, 2,4)=1. 

79 DCZN4=1. 

90 DDD=DCZNl=!'DCZN2’i‘DCZN3 
IF(OOD.NE.l .)G0 TO 20 
KCUNT=K0UNT+1 

IF(KOUNT.LE.NCOUNT)GO TO 10 
DC 100 1=1,4 

100 DATA( I, 1 ) = DATA( I , 1 )/NCOUNT 
CO 110 I=1,NC0UNT 

DC 110 11=1,4 

DATA( II,2) = (TEST( 1 , 1 , 1 1 ) -DATA ( I I , 1 ) ) ^#2 . /FLC AT ( NC GUNT- 
*1 )+DATA(II ,2) 

IF(T£ST(I,2,II).EC.l.) DATA ( 1 1 , 3 ) =DATA ( 11,3 )+l. 

110 CONTINUE 

DC 120 J=l,4 

120 DATA( J,4)=NC0UNT-DATA(J,3) 

WRITE(6,130) ((DATA( I , J ) , J = 1 , 4 ) , I =1 ,3) 

130 FORMAT!/ ,4X,‘ MEAN’ , FIO . 5 , • VAR’,F10.5,' NO ACCEPT • 

=5',F5.0,' NO REJECT ’,F5.0) 

101 FORMATdOFlO.S ) 

102 F0RMAT(5X,’ MEDIAN IS’,F7.3,’ SHAPE PARAMETER IS’, 
=i'F7.3) 

GO TO 300 
1000 STOP 
END 
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